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M.  ANaCraCI  (Umt  200  yarn »)  —  (  /'  ^ 

The  research  investigation  described  in  ^this  report  has  focused  on  a  preliminary  analy¬ 
sis  of  seismii  variability  at  the  Shagan  River  nuclear  test  site.  In  order  to  carry  out 
this  analysis ,m^arge  volumes  of  teleseismic  P  wave  amplitude  and  arrival  time  data  recorded 
from  explosions  at  this  test  site  were  collected  and  statistically  analyzed  in  an  attempt 
to  define  any  systematic  trends  Which  correlate  with  source  location.  The  results  of  the 
analysis  of  the  teleseismic  m£<data  have  been  shown  to  provide  strong  evidence  of  systema¬ 
tic  geophysical  variations  within  the  Shagan  River  testing  area.  In  particular,  the  mb 
station  corrections  have  been  shown  to  vary  with  source  location  in  a  manner  which  depends 
on  source-to-station  azimuth.  Results  of  detailed  comparisons  of  mg  residual  data  at  com¬ 
mon  stations  from  selected  pairs  of  nearby  explosions  have  been  used  to  conclude  that  the 
observed  variations  in  the  azimuthal  patterns  of  the  n(b)  residuals  are  related  to  event 
location  rather  than  tectonic  release  effects.  <5n  fact  th€  ^results  of  this  preliminary 
analysis  suggest  that  the  observed  m£  differences  may  be  associated  with  changes  in  the 
near-source  portions  of  the  teleseismic  P  wave  propagation  paths. 

In  contrast  to  the  residual  data,  the  reduced  P  wave  travel-time  residuals  were 
found  to  show  no  obvious  correlation  with  event  location,  and  this  has  been  interpreted  as 
an  indication  that  the  available  bulletin  data  are  not  precise  enough  to  resolve  the  rather 
small  travel-time  variations  which  might  be  expected  to  accompany  the  observed  mg  varia¬ 
tions. _ ;  •:  ■'  i<  ■"  ">■'  . _ _ 
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EXECUTIVE  SUMMARY 


This  report  summarizes  the  results  of  a  preliminary 
investigation  of  seismic  variability  within  the  Shagan  River 
region  of  the  Russian  nuclear  test  site  at  Semipalatinsk . 

The  primary  objective  of  this  study  has  been  to  assess  the 
feasibility  of  using  teleseismic  P  wave  amplitude  and  arri¬ 
val  time  data  recorded  from  explosions  at  this  test  site  to 
define  any  systematic  trends  which  correlate  with  source 
location . 

The  teleseismic  P  wave  data  base  which  has  been 
assembled  for  this  purpose  is  described  in  Section  II  where 
the  distributions  of  the  data  with  respect  to  source  and 
station  parameters  are  documented  in  detail.  These  data 
consist  of  individual  m^  readings  and  P  wave  arrival  times 
recorded  by  a  selected  network  of  94  worldwide  receiver  sta¬ 
tions  from  a  sample  of  52  Shagan  River  underground  nuclear 
explosions  which  have  been  assigned  m.  values  of  5.5  or 
greater.  These  seismic  data,  together  with  the  refined  ex¬ 
plosion  epicenters  and  origin  times  determined  by  Marshall 
et  al.  (1984) ,  have  provided  the  basis  for  investigating 
test  site  variability. 

Standard  seismological  and  statistical  procedures 
are  applied  to  these  data  in  Section  III  where  network- 
averaged  magnitudes  are  computed  for  each  explosion  and 
average  correction  factors  are  determined  for  each  of  the 
selected  stations  which  recorded  these  events .  The  varia¬ 
tions  of  these  m^  station  correction  factors  with  source 
location  are  then  carefully  analyzed  and  shown  to  strongly 
support  the  conclusion  that  systematic  geophysical  varia¬ 
tions  occur  within  the  Shagan  River  testing  area.  In  parti¬ 
cular,  it  is  demonstrated  that  the  corrected,  single  station 
residuals  for  explosions  in  the  southwest  portion  of  the 
test  site  appear  to  be  random  at  any  given  azimuth  and 


quite  similar  from  event  to  event,  while  there  are  large 
(0.5  units  m^)  variations  in  the  corrected  m^  residuals  with 
azimuth  between  explosions  in  close  proximity  in  the  north¬ 
east  and  central  portions  of  the  test  site.  Detailed  com¬ 
parisons  of  m^  residual  data  from  selected  pairs  of  explo¬ 
sions  are  then  used  to  demonstrate  that  the  observed  varia¬ 
tions  in  the  azimuthal  patterns  of  these  residuals  are  re¬ 
lated  to  event  location  rather  than  tectonic  release  effects. 
Moreover,  results  of  some  preliminary  analyses  are  presented 
which  suggest  that  the  variations  in  the  corrections  for 
stations  in  a  given  azimuth  are  systematic  enough  to  be  con¬ 
toured  as  a  function  of  source  location.  This  is  interpreted 
as  am  indication  that  the  observed  differences  are  associated 
with  changes  in  the  near-source  P  wave  propagation  paths  to 
teleseismic  distances  as  a  function  of  source  location  within 
the  test  site. 

The  results  of  some  preliminary  analyses  of  the  P 
wave  travel-time  data  are  also  presented  in  Section  III 
where  it  is  shown  that  the  reduced  P  wave  travel-time  resid¬ 
uals  do  not  show  the  correlation  with  event  location  inferred 
from  the  analysis  of  the  corresponding  residual  data. 

These  results  are  interpreted  as  an  indication  that  the  ISC 
travel  time  data  are  not  precise  enough  to  be  useful  for 
identifying  variations  in  the  subsurface  geology  across  the 
Shagan  River  test  site. 


I .  INTRODUCTION 


In  1974,  the  U.S.  and  U.S.S.R.  signed  a  Threshold 
Test  Ban  Treaty  (TTBT )  which  prohibits  the  testing  of  under¬ 
ground  nuclear  explosions  with  yields  greater  than  150  kilo- 
tons.  Upon  ratification,  the  treaty  calls  for  the  bilateral 
exchange  of  certain  geologic  and  geophysical  data,  as  well 
as  the  yields  of  two  calibration  events,  in  each  so-called 
"geophysically  distinct"  testing  area,  in  order  to  facilitate 
verification  of  treaty  compliance.  Although  not  defined 
explicitly  in  the  TTBT  protocol,  the  term  "geophysically 
distinct"  is  intended  to  denote  an  area  within  which  the 
geophysical  properties  controlling  the  magnitude-yield  rela¬ 
tionship  are  uniform  -  that  is,  an  area  within  which  a  sin¬ 
gle  yield-scaling  relation  holds  for  all  explosions.  Given 
adequate  calibration  data,  such  areas  are  generally  fairly 
easy  to  recognize.  For  example,  at  the  Nevada  Test  Site 
(NTS),  the  differences  between  the  seismic  coupling  charac¬ 
teristics  of  explosions  in  the  dry  alluvium  environment  of 
Yucca  Flat  versus  those  in  the  saturated  volcanic  environ¬ 
ment  of  Pahute  Mesa  are  easily  recognized  by  comparing  the 
magnitude-yield  populations  characteristic  of  those  two  test¬ 
ing  areas.  However,  for  areas  such  as  the  Semipalatinsk 
test  site  in  the  U.S.S.R.,  for  which  such  calibration  data 
are  not  available,  it  is  not  obvious  how  such  geophysically 
distinct  areas  could  be  identified  using  information  known 
to  us  at  the  present  time.  The  objective  of  the  preliminary 
research  investigation  described  in  this  report  has  been  to 
assess  the  feasibility  of  using  teleseismic  P  wave  data  re¬ 
corded  from  explosions  to  identify  geophysically  distinct 
testing  areas  within  the  Shagan  River  region  of  the  Semi¬ 
palatinsk  test  site. 

Figure  1  shows  the  geographic  location  of  the  prin¬ 
cipal  Russian  Nuclear  Testing  Ground  (NTG)  near  Semipalatinsk 
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Degrees  North  Latitude 


in  the  Eastern  Kazakhstan  Province  of  the  U.S.S.R.  It  can 

be  seen  that  the  NTG  is  quite  large,  covering  a  total  sur- 

2 

face  area  of  approximately  36,000  km  .  However,  most  of  the 

underground  nuclear  testing  conducted  at  this  site  occurs  in 

two  restricted  areas  known  as  the  Degelen  Mountain  and 

Shagan  River  testing  areas.  These  two  sites  are  outlined  in 

Figure  2  where  it  can  be  seen  that  the  Degelen  test  site  is 

located  in  the  central  part  of  the  NTG  and  occupies  a  sur- 

2 

face  area  of  approximately  410  km  while  the  larger  Shagan 

River  test  site  is  located  in  the  northeast  quadrant  of  the 

2 

NTG  and  occupies  a  surface  area  of  approximately  640  km  . 

In  recent  years,  the  largest  magnitude  Russian  underground 
nuclear  explosions  have  consistently  occurred  at  the  Shagan 
River  test  site  and,  consequently,  it  is  this  area  which  is 
currently  of  primary  interest  with  regard  to  yield  estima¬ 
tion.  The  contours  superimposed  on  the  maps  in  Figure  2 
denote  the  variations  in  crust  and  upper  mantle  structure 
across  the  area  as  inferred  from  Russian  gravity  and  magnetic 
data  (Zlavdinov,  1974)  .  It  can  be  seen  that  these  data  sug¬ 
gest  some  systematic  variations  in  crustal  structure  across 
the  NTG.  Thus,  for  example,  in  the  vicinity  of  the  Degelen 
site  the  sedimentary  cover  is  virtually  absent,  the  granitic 
.tayer  has  an  average  thickness  of  20  km,  the  basaltic  layer 
has  a  maximum  thickness  of  21  km  and  the  depth  to  the  Moho 
is  as  great  as  41  km;  while  at  the  Shagan  River  site  the 
sedimentary  layer  thickens  to  as  much  as  7  km,  the  granitic 
and  basaltic  layers  both  decrease  in  thickness  (to  13  km 
and  18.5  km  respectively)  and  the  average  Moho  depth  is 
about  39  km.  It  is  reasonable  to  expect  that  these  large 
scale  crustal  variations  may  be  associated  with  smaller 
scale  variations  in  structure  at  the  shallower  depths  where 
testing  is  taking  place,  and  it  is  these  latter  variations 
which  we  will  be  attempting  to  map  through  comparisons  of 
the  teleseismic  P  waves  recorded  from  various  Shagan  River 
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NE-SW  hatched  pattern  and  Degelen  Mtn.  (central) 
hatched  pattern.  (Adapted  from  Zlavdinov,  1974) . 


underground  explosions. 

The  organization  of  this  report  may  be  briefly  sum¬ 
marized  as  follows.  The  data  samples  of  teleseismic  P  wave 
amplitudes  and  travel  times  compiled  for  the  purposes  of 
this  study  are  described  in  detail  in  Section  II.  In  Sec¬ 
tion  III  these  data  are  analyzed  using  a  variety  of  statis¬ 
tical  analysis  techniques  in  an  attempt  to  isolate  any  signi 
ficant  correlations  with  event  location.  This  is  followed 
in  Section  IV  by  a  summary  and  a  listing  of  preliminary  con¬ 
clusions  concerning  variations  in  teleseismic  P  wave  charac¬ 
teristics  with  source  location  within  the  Shagan  River  test 


II.  DESCRIPTION  OF  THE  TELESEISMIC 
P  WAVE  DATA  BASE 

The  data  sample  which  has  been  compiled  to  investi¬ 
gate  variability  at  the  Shagan  River  test  site  is  composed 
of  teleseismic  P  wave  observations  recorded  from  underground 
nuclear  explosions  conducted  at  that  site.  In  assembling 
this  sample,  only  those  explosions  which  occurred  in  the 
time  interval  1964  through  1982  which  have  been  assigned 
values  of  greater  than  5.5  by  ISC  have  been  considered. 

This  time  interval  was  selected  to  coincide  with  that  for 
which  data  are  available  from  either  the  ISC  Bulletins  or  the 
NEIS  Earthquake  Data  Reports,  while  the  m^  threshold  of  5.5 
was  adopted  to  guarantee  reasonably  uniform  recording  by  the 
worldwide  networks  of  stations.  Using  these  criteria,  a 
sample  of  52  Shagan  River  explosions  has  been  identified  for 
analysis.  The  ISC  source  parameters  for  these  52  explosions 
are  listed  in  Table  1,  together  with  the  epicenter  locations 
obtained  for  these  same  events  by  Marshall  et  al.  (1984) 
using  the  Joint  Epicenter  Determination  (JED)  method.  These 
latter  locations  are  believed  to  be  considerably  more  accu¬ 
rate  in  that  they  have  been  computed  relative  to  the  known 
location  of  the  Shagan  River  cratering  explosion  of  1/15/65 
(event  #  1  in  Table  1)  and,  consequently,  will  be  used  for 
the  purposes  of  this  investigation.  The  JED  locations  of 
the  52  explosions  of  Table  1  are  displayed  graphically  in 
Figure  3  where  it  can  be  seen  that  they  are  fairly  broadly 
distributed  across  the  Shagan  River  testing  area.  The  95 
percent  confidence  ellipses  about  these  epicenters  reported 
by  Marshall  et  al .  (1984)  suggest  that  these  locations  are 
accurate  to  within  a  few  kilometers  in  most  cases.  The  dis¬ 
tribution  of  these  events  with  respect  to  ISC  value  is 
shown  in  Figure  4  where  it  can  be  seen  that  the  majority  of 
the  explosions  have  been  assigned  values  in  the  range 
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5.6  £  i  6.2,  which  is  well  within  the  optimum  operating 
range  of  the  worldwide  seismic  networks  which  report  to  ISC. 
At  the  present  time,  the  source  parameters  for  all  events 
occurring  since  8/31/82  are  available  from  the  NEIS  service 
only;  NEIS  m^  values  differ,  on  average,  by  0.10  m^  units  or 
less  from  ISC  m^  values  for  this  source  region. 

He  have  extracted  all  available  teleseismic  P  wave 
observations  for  each  selected  explosion  from  earthquake 
catalogs  published  by  the  ISC  and  NEIS  services.  Over  the 
time  period  encompassing  our  explosion  sample  (i.e.,  1964 
through  1982) ,  a  total  of  approximately  863  different  world¬ 
wide  seismic  stations  and  arrays,  ranging  in  epicentral  dis¬ 
tance  from  30°  to  105°  and  widely  distributed  in  azimuth, 
reported  teleseismic  P  wave  detections  from  some  of  the  se¬ 
lected  52  Shagan  events.  Of  course,  many  of  the  863  world¬ 
wide  receiver  stations  were  operational  only  within  specific 
time  periods  and,  thus,  the  reporting  consistency  varies 
significantly  between  stations.  Our  aim  in  defining  an 
optimum  data  base  has  been  to  preserve  the  maximum  number  of 
common  teleseismic  stations,  with  both  wide  azimuthal  and 
epicentral  distance  coverage,  having  detected  all  or  most  of 
the  selected  Shagan  River  explosions  and  reported  all  avail¬ 
able  event  parameters .  The  event  parameters  that  are  most 
useful  for  our  statistical  analyses  are  P  wave  arrival  time, 
travel  time  residual  computed  with  respect  to  the  travel  time 
expected  for  a  given  distance  in  a  standard,  average  earth 
model  and  individual  station  m^,  or  the  log  of  the  ratio  of 
P  wave  amplitude  to  period  (i.e.,  log  A/T) .  Every  station 
detecting  a  given  explosion  in  the  event  sample  has  reported 
at  least  P  wave  arrival  time  information  allowing  for  the 
determination  of  travel  time?  therefore,  the  detection  cri¬ 
teria  for  the  selection  of  a  suitable  network  of  stations 
used  in  this  study  has  been  based  primarily  on  observed 
arrival  time  data.  As  such,  we  have  found  that  the 


meaningful  station  coverage  consists  of  94  worldwide  stations 
having  detected  at  least  75  percent  of  the  total  explosion 
sample.  The  mean  detection  level  for  this  group  of  stations 
is,  however,  87  percent.  There  appears  to  be  no  significant 
advantage  to  selecting  a  detection  level  higher  or  lower 
than  75  percent  since  there  is  at  least  a  25  percent  decrease 
in  the  station  sample  if  the  event  detection  level  is  con¬ 
strained  to  80  percent  or  higher;  and,  although  a  large  in¬ 
crease  (38  percent)  in  the  total  station  sample  is  achieved 
by  reducing  the  detection  requirement  to  65  percent,  there 
is  no  net  gain  in  the  range  of  recording  azimuths. 

Unfortunately,  due  to  inconsistency  in  station  re¬ 
porting,  the  number  of  stations  for  which  individual  station 
data  are  available  is  frequently  significantly  smaller 
than  the  total  number  reporting  arrival  time  data.  By  merg¬ 
ing  station  data  published  in  both  ISC  and  NEIS  catalogs,  it 
has  been  possible  to  extract  the  maximum  number  of  station 
(or  log  A/T)  values.  Of  the  final  94  stations,  78  (i.e., 

83  percent)  have  reported  ( or  log  A/T)  values  for  the 

Shagan  explosion  sample.  The  mean  number  of  m^  (or  log  A/T) 
data  available  for  this  group  is  27  ±  3  out  of  a  possible 
52  events,  or  52  percent  of  the  entire  explosion  sample. 

The  statistical  significance  of  the  size  of  this  m^  data  base, 
relative  to  the  size  of  the  travel  time  data  base,  will  be 
addressed  in  Section  III.  In  cases  where  log  A/T  quantities 
have  been  reported,  m^  values  were  calculated  based  on  the 
appropriate  Gutenberg-Richter  distance  correction  factors. 

As  a  result  of  the  above  descriptions  of  the  complete  data 
base  utilized  in  this  study,  we  find  that  for  our  selected 
seismic  network  of  94  stations  there  are,  on  the  average, 

82  stations  (or  a  range  of  31  to  94)  with  travel  time  infor¬ 
mation  and  41  stations  (or  a  range  of  15  to  59)  with  individ¬ 
ual  station  m^  (or  log  A/T)  values  reported  per  event.  To 
ensure  effective  statistical  and  seismological  analyses,  we 


have  compiled  these  data  consisting  of  several  thousand  in¬ 
dividual  station  arrival  time  and  m^  values,  carefully  veri 
fied  and  assembled  into  a  tractable  data  base  on  the  DARPA 
Center  for  Seismic  Studies  (CSS)  computer  system  in  Rosslyn 
Virginia. 

The  individual  station  information  that  is  perti¬ 
nent  to  this  investigation,  including  epicentral  distance 
and  recording  azimuth  to  the  center  of  the  Shagan  test  site 
as  well  as  number  of  event  recordings  available,  are  listed 
in  Table  2.  Figure  5  shows  the  azimuthal  distribution  of 
the  teleseismic  station  coverage  common  to  the  Shagan  explo 
sion  sample  where  it  can  be  seen  that,  although  there  is 
generally  wide  azimuthal  coverage,  the  predominant  coverage 
occurs  in  the  west-northwest  and  north -northeast  azimuthal 
sectors.  The  distribution  of  these  same  stations  on  the 
basis  of  epicentral  distance  to  the  Shagan  test  site  is  dis 
played  in  Figure  6  indicating  that  the  majority  of  stations 
are  between  30  and  60  degrees  away  from  Shagan.  The  actual 
map  locations  of  the  stations  are  shown  on  an  azimuthally 
equidistant  projection  centered  on  the  Shagan  River  test 
site  in  Figure  7,  where  it  can  be  seen  that  the  heaviest 
concentration  of  stations  occurs  in  western  Europe. 


Table  2 


Station 

Code 

A° 

TELESEISMIC 

STATION  INFORMATION 

Number  of  Event  Detections 

AZ° 

n^Cor  log  A/T) 

P-Wave 

Arrival  Time 

ADK 

60.1 

44 

5 

44 

AKU 

48.6 

327 

28 

46 

ALE 

46.4 

353 

35 

41 

ALQ 

95.4 

4 

31 

52 

ASP 

88.0 

131 

21 

44 

BER 

41.1 

314 

1 

47 

BKS 

90.6 

17 

36 

48 

BRG 

40.2 

297 

45 

48 

BSE 

45.8 

297 

32 

48 

BUD 

38.7 

290 

32 

40 

BUL 

82.7 

227 

32 

51 

CDF 

45.2 

297 

32 

50 

CHG 

35.0 

145 

24 

41 

CLL 

40.5 

298 

48 

52 

COL 

59.9 

21 

25 

45 

COP 

39.2 

305 

44 

47 

CTA 

91.5 

120 

25 

47 

DAG 

43.8 

341 

37 

45 

DOU 

45.8 

300 

5 

52 

EDM 

76.7 

8 

31 

42 

ERA 

47.4 

310 

33 

48 

EUR 

90.1 

11 

26 

43 

FFC 

75.8 

1 

31 

41 

FHC 

87.4 

17 

0 

41 

fln 

49.4 

301 

24 

47 

FRI 

91.9 

15 

0 

48 

FUR 

43.0 

295 

34 

50 

15 


Table  2  (Cont'd) 


Number  of  Event  Detections 


Station  P-Wave 

Code  A°  AZ°  m^or  lo9  A/T)  Arrival  Time 


10  Stations 


20  Stations 


30  Stations 


40  Stations 
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Azimuthal  station  coverage  for  the  initial  station 
network  used  in  this  study. 
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Figure  6.  Histogram  showing  the  frequency  distribution  of 
network  stations  by  epicentral  distance  to  the 
center  of  the  Shagan  test  site. 
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Figure  7.  Map  locations  of  stations  shown  on  an  azimuthally 
equidistant  projection  centered  on  the  Shagan 
River  test  site,  denoted  by  solid  triangle  (A)  . 


III.  ANALYSIS  OF  THE  TELESEISMIC  P  HAVE  DATA 


3.1  DATA 

The  m^  data  base  which  was  described  in  detail  in 
the  preceding  section  has  been  subjected  to  a  variety  of 
statistical  analyses  in  an  attempt  to  identify  any  systema¬ 
tic  geophysical  variations  within  the  Shagan  River  testing 
area.  For  the  purposes  of  this  investigation,  the  individ¬ 
ual  station  m^  values  have  been  initially  processed  using 
the  LSMF  (least  squares  matrix  factorization)  procedure  orig¬ 
inally  described  by  Douglas  (1966) .  Under  this  formulation, 
data  recorded  on  a  network  of  stations  from  explosions  at  a 
given  test  site  are  statistically  analyzed  to  simultaneously 
determine  the  best-fit  station  correction  factors  and  net¬ 
work-averaged  magnitudes  for  each  event,  under  the  constraint 
that  the  station  correction  factors  stun  to  zero.  The  prin¬ 
cipal  advantage  of  this  approach  is  that  it  provides  a  con¬ 
sistent  framework  for  analyzing  data  recorded  on  networks 
of  stations  for  which  different  stations  have  recorded  dif¬ 
ferent  numbers  of  the  events  under  consideration.  In  the 
context  of  the  present  study,  it  is  important  to  note  that 
the  "station  corrections"  computed  using  the  LSMF  procedure 
actually  represent  the  composite  of  effects  at  the  source, 
along  the  propagation  path  and  at  the  receiver  which  cause 
the  values  at  a  particular  station  to  be  consistently 
different  from  the  corresponding  large  network  average 
values.  Thus,  for  explosions  at  a  particular  test  site, 
the  "station  corrections"  may  be  more  closely  associated 
with  propagation  effects  near  the  source  than  with  varia¬ 
tions  in  the  crustal  structure  beneath  the  receivers.  A 
dramatic  example  of  this  fact  is  provided  in  Figure  8  which 
shows  a  comparison  of  network-averaged  values  for  selected 
explosions  at  the  Degelen  Mountain  and  Shagan  River  testing 
areas  of  the  Semipalatinsk  Test  Site  with  the  corresponding 
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Figure  8.  Comparison  of  single  station  HFS  and  PDE  network 
mb  values  for  underground  nuclear  explosions  at 
the  Degelen  and  Shagan  River  test  sites. 


m^  values  observed  at  station  HFS  (Hagfors)  in  Sweden.  Note 
that  the  data  for  the  two  areas  segregate  into  two  dis¬ 
joint  populations  such  that,  for  a  fixed  value  of  network 
m^,  the  HFS  m^  values  are  nearly  a  full  order  of  magnitude 
larger  for  explosions  at  Shagan  them  for  explosions  at 
Degelen.  That  is,  the  HFS  m^  "station  correction”  factor 
with  respect  to  the  worldwide  average  changes  by  nearly  a 
full  order  of  magnitude  between  these  two  test  sites  which 
are  separated  by  less  than  60  km,  presumably  reflecting  near¬ 
source  propagation  path  effects  associated  with  the  well- 
known  differences  in  the  near-surface  geologic  environments 
at  the  two  sites  (i.e.,  granitic  intrusives  at  Degelen  versus 
sedimentary  formations  at  Shagan) .  This  example  demonstrates 
that  teleseismic  m^  variations  do  correlate  with  variations 
in  test  site  conditions,  at  least  in  some  cases. 

The  analysis  of  the  Shagan  River  data  set  was  ini¬ 
tiated  by  running  2m  LSMF  analysis  on  the  individual  station 
m^  data  recorded  from  the  selected  52  explosions  described 
in  Section  II.  A  total  of  71  stations  which  reported 
values  for  5  or  more  of  these  explosions  were  included  in 
this  preliminary  evaluation.  This  number  was  subsequently 
reduced  to  59  by  eliminating  those  stations  for  which  the 
computed  standard  error  of  estimate  in  the  mean  station  cor¬ 
rection  exceeded  0.3  magnitude  units.  This  judgment  was 
made  on  the  basis  of  the  observation  that  the  average  stand¬ 
ard  error  of  estimate  in  the  mean  station  correction  for  the 
final  59  stations  is  less  than  0.2  magnitude  units.  The 
number  of  events  for  which  values  were  reported,  for  each 
of  these  59  stations, are  displayed  in  histogram  form  in  Fig¬ 
ure  9  where  it  can  be  seen  that  all  but  7  of  these  stations 
reported  values  for  at  least  20  of  the  52  selected  Shagan 
River  explosions.  An  LSMF  analysis  was  run  on  this  final 
data  set  to  obtain  the  network-averaged  m^  values  and  aver¬ 
age  station  corrections  to  be  used  in  the  subsequent 
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investigations.  The  results  of  these  calculations  were 
found  to  be  quite  consistent  with  those  reported  by  Marshall 
et  al.  (1984)  on  the  basis  of  their  analysis  of  data  recorded 
from  a  similar  set  of  Shagan  River  explosions  on  a  large, 
worldwide  network  consisting  of  174  stations,  of  which  54  are 
in  common  with  the  network  used  in  the  present  study.  This 
consistency  is  graphically  illustrated  in  Figures  10  and  11 
where  the  network-averaged  m^  values  and  station  corrections 
estimated  in  the  present  study  are  compared  with  those  of 
Marshall  et  al.  (1984)  .  Complete  lists  of  the  final  network 
m ^  values  and  station  corrections  derived  in  this  study  are 
presented  in  Tables  3  and  4  respectively.  For  purposes  of 
comparison,  the  corresponding  network  m^  estimates  obtained 
by  Marshall  et  al.  (1984)  and  the  ISC  are  also  included  in 
Table  3  for  those  explosions  for  which  they  are  available. 

Having  obtained  average  station  corrections  for 
explosions  at  the  Shagan  River  site,  it  now  remains  to 
examine  variability  within  the  test  site  itself.  We  have 
attempted  to  accomplish  this  through  systematic  evaluation 
of  the  observed  distributions  of  station-corrected  m^  resid¬ 
uals  (i.e.,  corrected  single  station  m^  -  event  n^)  as  a 
function  of  source  location.  Now,  in  the  absence  of  any 
underlying  deterministic  mechanisms,  it  would  be  expected 
that  the  residuals  for  a  given  station  would  be  randomly 
distributed  with  respect  to  event  location.  However,  this 
is  not  found  to  be  the  case  at  Shagan  River.  For  example. 
Figure  12  shows  the  station-corrected  m^  residuals  as  a 
function  of  event  location  for  stations  MAT  and  ERA.  Note 
that  the  residuals,  which  can  be  regarded  as  variations  in 
the  "station  corrections"  with  source  location,  show  pro¬ 
nounced  trends  in  that  residuals  of  the  same  size  and  sign 
tend  to  cluster  into  geographical  groups.  This  suggests 
that  there  are  some  unaccounted  for  source  region  physical 
mechanisms  which  are  affecting  the  radiation  of  P  wave 


m £  Station  Correction  (This  Study) 


Comparison  of  station  corrections  derived  in  this 
study  with  those  of  Marshall  et  al.  (1984)  for  com¬ 
mon  stations. 
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Table  4 


STATION  CORRECTIONS 


Station 

This  Study 
(#  of  Events) 

Marshall  et  al 
(#  of  Events) 

ALE 

-0.404 

(35) 

-0.416 

(32) 

ALQ 

-0.483 

(31) 

- 

BKS 

0.075 

(36) 

0.077 

(33) 

BRG 

-0.458 

(45) 

-0.533 

(40) 

BSF 

-0.179 

(32) 

-0.159 

(29) 

BUL 

0.018 

(32) 

0.060 

(44) 

CDF 

-0.618 

(32) 

-0.604 

(31) 

CHG 

-0.111 

(24) 

-0.005 

(20) 

CLL 

0.016 

(48) 

-0.061 

(42) 

COL 

0.399 

(25) 

0.315 

(22) 

COP 

-0.008 

(44) 

-0.019 

(36) 

CTA 

-0.132 

(25) 

-0.141 

(7) 

DAG 

-0.143 

(37) 

-0.280 

(37) 

DOU 

0.372 

(5) 

- 

EDM 

0.377 

(31) 

0.371 

(30) 

ERA 

0.170 

(33) 

0 .153 

(39) 

FFC 

0.325 

(31) 

0.347 

(33) 

FLN 

0.289 

(24) 

0.343 

<21 ) 

GDH 

0  .016 

(22) 

0.031 

(31) 

GRF 

-0.118 

(35) 

-0.112 

(33) 

GRR 

-0.117 

(31) 

-0.097 

(25) 

HAU 

0.003 

(34) 

0.037 

(32) 

HES 

0.993 

(34) 

0.937 

(43) 

HOF 

-0  .057 

(19) 

-0.106 

(14) 

HYB 

0.120 

(45) 

0.125 

(37) 

INK 

0.501 

(28) 

0.393 

(26) 

KHC 

-0.272 

(41) 

-0.329 

(37) 
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KIR 

0.970 

(12) 

— 

KJF 

0.582 

(44) 

0.471 

(43) 

KRA 

0.242 

(43) 

0.122 

(40) 

LBF 

-0.310 

(28) 

-0.293 

(26) 

LFF 

-0.132 

(28) 

-0.101 

(26) 

LJU 

-0.096 

(7) 

- 

LMR 

-0.182 

(30) 

-0.071 

(27) 

LOR 

-0.096 

(30) 

-0.071 

(34) 

LPF 

0.025 

(24) 

0.017 

(20) 

LPO 

0.145 

(25) 

0.149 

(21) 

LRG 

-0.087 

(34) 

-0.029 

(30) 

LSF 

-0.466 

(29) 

-0.501 

(26) 

MAT 

-0.711 

(30) 

-0.762 

(16) 

MFF 

0.045 

(27) 

0.014 

(24) 

MNT 

-0.132 

(24) 

-0.120 

(18) 

MOX 

-0.072 

(48) 

-0.072 

(41) 

MUN 

0.303 

(13) 

0.290 

(ID 

NEW 

-0.459 

(15) 

-0.481 

(14) 

OTT 

0.123 

(33) 

0.167 

(30) 

PMG 

0.254 

(23) 

0.287 

(22) 

PNT 

-0.172 

(33) 

-0.211 

(32) 

PRU 

-0.513 

(43) 

-0  .533 

(40) 

RJF 

-0.317 

(30) 

-0.331 

(27) 

SES 

0.391 

(26) 

0.312 

(27) 

SSC 

0.097 

(29) 

0.183 

(28) 

SSF 

-0.326 

(31) 

-0.274 

(29) 

STU 

0.003 

(31) 

-0.062 

(32) 
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Figure  12.  Station-corrected  nik  residuals  as  a  function  of  Shagan  event  location 
for  stations  MAT  and  EKA  at  the  denoted  recording  azimuths. 


energy  to  teleseismic  distances.  In  an  attempt  to  examine 
these  effects  in  more  detail,  we  have  analyzed  the  azimuthal 
distribution  of  the  corrected  m^  residuals  on  an  event-by¬ 
event  basis.  Figure  13  shows  the  results  for  four  typical 
southwest  Shagan  events.*  In  these  and  subsequent  polar 
coordinate  plots  (station  azimuth  measured  clockwise  from 
north) ,  the  circle  corresponds  to  a  corrected  residual  of 
zero,  while  positive  residuals  plot  outside  the  circle  and 
negative  residuals  plot  inside  the  circle  according  to  the 
-0.5  to  +0.5  m^  unit  scale  shown  on  the  figure.  Note  from 
Figure  13  that  for  these  explosions  in  the  southwest  portion 
of  the  test  site,  the  corrected  residuals  appear  to  be 
random  at  any  given  azimuth,  and  quite  similar  from  event- 
to-event.  However,  as  is  illustrated  in  Figure  14,  this 
simple  picture  does  not  hold  for  explosions  in  the  northeast 
and  central  portions  of  the  test  site.  It  can  be  seen  that 
in  these  areas  there  are  some  large  variations  (i.e.,  ±  0.5 
units)  in  the  corrected  m^  residuals  with  azimuth  between 
events  in  very  close  proximity.  On  a  more  detailed  level, 
it  is  difficult  to  quantitatively  assess  the  magnitude  of 
these  variations  from  comparisons  such  as  those  shown  in 
Figure  14  in  that  somewhat  different  networks  of  stations 
recorded  each  of  these  events.  Therefore,  in  an  attempt  to 
eliminate  any  network  biasing  effects,  we  have  also  compared 
corrected  m^  residuals  for  different  pairs  of  events  at 
common  sets  of  stations .  Four  such  comparisons  are  shown  in 
Figures  15-18  for  selected  explosions  located  in  the  same 
area  covered  by  Figure  14.  It  can  be  seen  that  these  in¬ 
dividual  event  comparisons  at  common  stations  are  quite 
consistent  with  the  general  trends  shown  in  Figure  14,  con¬ 
firming  the  fact  that  large  variations  in  corrected  m^ 

*  The  event  numbers  (#)  shown  on  this  and  subsequent  figures 
are  referenced  to  the  event  sequence  numbers  given  in  Table 
3. 
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Figure  14.  Azimuthal  distribution  of  corrected  m^  residuals  as  a  function  of 
event  location;  NE  and  central  events.  Circle  represents  zero 
residual;  range  -.5  m.  to  +.5  m.  is  denoted. 


Figure  15.  Comparison  of  corrected  mb  residuals  for  events  #52  and  #14  at 
34  common  stations. 


3  common  stations. 


residuals  do  occur  between  events  in  close  proximity  in  the 
northeast  and  central  portions  of  the  Shagan  River  test  site. 
Specifically,  explosions  in  the  northeast  quadrant  (e.g., 
events  #52,  #25,  and  #6)  show  large  positive  corrected  m^ 
residuals  at  stations  in  the  northern  azimuth  window  and 
predominantly  negative  corrected  m^  residuals  at  stations  in 
the  northwest  azimuth  window,  while  explosions  located  nearer 
to  the  center  of  the  test  site  (e.g.,  events  #14  and  #30) 
show  opposite  trends,  with  large  negative  corrected  m^  resid¬ 
uals  at  stations  in  the  northern  azimuth  and  predominantly 
positive  corrected  residuals  at  stations  in  the  northwest 
azimuth.  Moreover,  these  differences  are  large  enough  that 
they  are  statistically  significant  at  a  very  high  confidence 
level.  This  fact  is  illustrated  in  Figure  19  which  shows 
the  differences  in  the  corrected  m^  residuals  at  common  sta¬ 
tions  between  events  #52  and  #14  plotted  as  a  function  of 
azimuth.  Note  that  the  average  m-,  residual  difference  com¬ 
puted  for  stations  in  the  260°  to  300°  azimuth  window  is 
significantly  different  from  that  computed  for  stations  in 
the  340°  to  20°  azimuth  window,  amounting  to  nearly  0.6  magni¬ 
tude  units.  It  follows  that  the  network-averaged  values 
and  corresponding  yield  estimates  determined  for  these  ex¬ 
plosions  will  be  highly  dependent  on  the  specific  azimuthal 
distributions  of  the  stations  used  to  compute  the  averages. 
Thus,  since  yield  estimation  is  generally  accomplished  using 
fixed  networks  of  stations,  some  consideration  must  be  given 
to  the  possible  effects  of  network  bias  on  the  yield  estimates 
for  explosions  in  this  portion  of  the  Shagain  River  test  site. 

Given  that  the  station-corrected  teleseismic 
data  provide  strong  evidence  of  systematic  geophysical  vari¬ 
ations  within  the  Shagan  River  testing  area,  it  now  remains 
to  attempt  to  identify  specific  causes  of  these  variations. 

In  the  light  of  the  evidence  presented  above ,  there  are  two 
alternate  hypotheses  which  seem  worthy  of  consideration: 
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Figure  19.  Station-corrected  residual  differences  between  event 
#52  and  event  #14  as  a  function  of  station  azimuth. 
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(1)  that  tectonic  release  induced  by  the  explosions  is  af¬ 
fecting  the  short-period  P  waves  and  (2)  that  there  are  var¬ 
iations  in  the  near-source  P  wave  propagation  paths  to  tele- 
seismic  distances  as  a  function  of  source  location  within 
the  test  site.  In  principle,  it  should  be  fairly  easy  to 
discriminate  between  these  two  very  different  hypotheses. 

In  fact,  however,  due  to  limitations  in  the  distribution  of 
teleseismic  P  wave  observations  with  azimuth  and  epicentral 
distance,  it  is  not  easy  to  definitively  test  either  hypoth¬ 
esis.  For  this  reason,  we  have  selected  a  control  group  of 
four  nearby  explosions  in  the  northeast  sector  of  the  Shagan 
River  testing  area  to  be  investigated  in  detail.  The  rela¬ 
tive  locations  of  these  four  explosions  (e.g.,  #25,  #28, 

#15  and  #41,  Cf.  Table  3)  within  the  Shagan  River  testing 
area  are  shown  in  Figure  20  where  it  can  be  seen  that  no  two 
of  the  events  are  separated  by  more  than  5  km.  The  letters 
A,  B  and  C  in  parentheses  beside  the  event  numbers  on  this 
figure  denote  the  tectonic  release  classification  of  these 
events  assigned  by  North  and  Fitch  (1981)  on  the  basis  of 
analyses  of  surface  wave  recordings  from  these  events ,  where 
A  denotes  an  explosion  accompanied  by  a  low  level  of  long 
period  Love  wave  excitation  and  no  observed  Rayleigh  wave 
phase  reversals,  B  denotes  an  explosion  with  observed  Ray¬ 
leigh  wave  phase  reversals  in  some  azimuths  and  C  denotes 
an  explosion  with  observed  Rayleigh  wave  phase  reversals  at 
most  stations.  That  is,  in  terms  of  the  long-period  surface 
waves  an  "A-type"  explosion  is  interpreted  to  be  one  accom¬ 
panied  by  a  low  level  of  tectonic  release,  while  a  "C-type" 
explosion  shows  convincing  evidence  of  very  strong  tectonic 
release  effects.  Thus,  as  can  be  seen  from  Figure  20,  the 
selected  four  events  are  not  only  closely  spaced  within  the 
region  where  the  azimuthal  patterns  of  the  corrected  m^ 
residuals  are  changing  rapidly,  but  also  represent  all  three 
surface  wave  tectonic  release  classifications.  it  follows 
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under  consideration. 

Regardless  of  which  of  the  two  hypotheses  is  cor¬ 
rect,  there  should  be  some  distance  dependence  to  the  cor¬ 
rected  m^  residual  pattern  at  a  fixed  azimuth  associated 
with  variations  in  P  wave  take-off  angle  at  the  source.  A1 
though  the  available  data  are  not  ideally  distributed  to  de 
fine  the  details  of  such  a  distance  dependence,  they  do  see 
to  confirm  its  existence.  For  example.  Figure  21  shows  the 
corrected  residuals  for  the  four  events  of  Figure  20, 
plotted  as  a  function  of  epicentral  distance.  A,  using  the 
data  recorded  at  the  30  stations  located  in  the  narrow  azi¬ 
muth  window  290°  <  0  <  305°.  In  both  this  and  the  next  few 
figures,  the  stations  are  shown  evenly  spaced  in  order  of 
increasing  azimuth  along  the  abscissa  to  permit  the  data 
from  these  tightly  clustered  stations  to  be  differentiated 
from  one  smother  (i.e.,  the  azimuth  scale  is  nonlinear). 

The  signs  and  sizes  of  the  observed  residuals  are  shown  in 
terms  of  squares  (positive)  and  triangles  (negative)  of 
varying  sizes,  keyed  to  the  scale  shown  at  the  top  of  the 
figure.  It  can  be  seen  from  these  four  examples  that  there 
is  indeed  a  pronounced  distance  dependence  to  the  corrected 
m^  residuals  at  this  northwest  azimuth,  characterized  by  a 
consistent  change  in  sign  of  the  residuals  at  an  epicentral 
distance  of  about  45  degrees.  Another  noteworthy  point  is 
that  the  residual  patterns  seem  to  correlate  better  with 
event  location  than  with  tectonic  classification.  That  is, 
the  patterns  for  the  two  northern  events  (#25  and  #28)  show 
at  the  top  of  Figure  21  are  quite  similar  despite  the  dif¬ 
ference  in  tectonic  classification  (i.e.,  A  versus  C)  and 
different  from  those  associated  with  the  two  southern  event 
shown  at  the  bottom  of  the  figure  which  have  been  assigned 
tectonic  classifications  of  B  (#41)  and  C  (#15)  . 


Figure  21.  Corrected  n^,  residuals  for  events  in  Figure  20  as  a  function  of  epicentral 
distance  determined  for  stations  within  the  azimuthal  range  290°  <  8  <  305 


Specifically,  the  signs  of  the  residuals  at  individual  sta¬ 
tions  change  in  going  from  north  to  south,  but  not  with 
tectonic  classification.  Moreover,  it  is  also  important  to 
note  that,  for  each  of  these  four  events,  the  average  resid¬ 
ual  in  the  epicentral  distance  range  of  45°  <  A  <  50°  differs 
from  that  in  the  epicentral  distance  range  40°  <  A  <  45°  by 
about  0.3  to  0.4  magnitude  units,  an  amount  which  is  signi¬ 
ficantly  greater  them  the  approximately  0.1  magnitude  shift 
that  would  be  theoretically  predicted  over  this  distance 
range  for  the  generally  inferred  mode  of  tectonic  release 
at  Shagan  (i.e.,  thrust  motion  on  a  fault  plane  dipping  at 
about  45°)  .  Although  these  results  argue  against  a  tectonic 
release  explanation  for  the  observed  m^  variability,  they 
are  not  really  definitive  in  that  they  are  representative 
of  only  a  single  narrow  azimuth  window.  Unfortunately,  the 
only  other  azimuth  window  for  which  there  is  enough  data  to 
make  meaningful  comparisons  is  the  northern  one  which  covers 
the  range  340*  <  6  <  20°.  The  corrected  residuals  for 
the  four  events  of  Figure  20  are  plotted  as  a  function  of 
epicentral  distance  in  Figure  22  using  the  data  recorded  at 
the  13  stations  located  in  this  northern  azimuth  window.  It 
can  be  seen  that  at  this  azimuth  data  are  available  only 
over  the  epicentral  distance  range  45*  <  A  <  95°  and  that, 
over  this  range,  there  is  no  clear  distance  dependence  to 
the  corrected  m^  residuals.  However,  as  with  the  north¬ 
western  azimuth  results  shown  in  Figure  21,  these  residual 
patterns  again  appear  to  correlate  better  with  event  loca¬ 
tion  than  with  tectonic  classification,  with  the  two  north¬ 
ern  events  showing  predominantly  positive  residuals  and  the 
two  southern  events  showing  predominantly  negative  residuals 
over  this  distance  range. 

Perhaps  a  more  direct  way  of  comparing  the 
data  for  the  four  events  of  Figure  20  is  to  compute  the 
differences  between  the  m^  residuals  at  common  stations. 
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Figure  22.  Corrected  m^  residuals  for  events  in  Figure  20 
as  a  function  of  epicentral  distance  determined 
for  stations  within  the  azimuthal  range  340°  <  9  < 
20®. 
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first  for  events  with  similar  tectonic  classification  and 
then  for  events  with  similar  locations.  Figure  23  shows  the 
differences  in  the  m^  residuals  at  common  stations  plotted 
as  a  function  of  station  azimuth  for  the  two  sets  of  events 
with  similar  tectonic  classification  (i.e.,  #15  versus  #28 
and  #41  versus  #25)  .  Note  that  in  both  these  cases  there  is 
evidence  of  pronounced  azimuthal  differences  between  the 
pairs  of  events  with  large  positive  residual  differences  to 
the  northwest  and  large  negative  residual  differences  to  the 
north.  In  contrast  to  this.  Figure  24  shows  the  m^  residual 
differences  at  common  stations  computed  for  the  pairs  of 
northern  (i.e.,  #28  versus  #25)  and  southern  (i.e.,  #41 
versus  #15)  events.  In  this  case,  the  residual  differences 
shows  no  pronounced  azimuthal  dependence  and  are  smaller, 
despite  the  fact  that  event  #28  (C)  is  at  a  greater  distance 
from  event  #25  (A)  than  from  event  #15  (C) .  These  results 
strongly  suggest  that  the  variations  in  the  azimuthal  pat¬ 
terns  shown  in  Figures  13-18  are  related  to  event  location 
rather  them  tectonic  release  effects  and,  thus,  support  the 
hypothesis  that  they  are  associated  with  variations  in  the 
near-source  P  wave  propagation  paths  to  teleseismic  dis¬ 
tances  . 

Considering  the  test  site  as  a  whole,  if  the  ob¬ 
served  azimuthal  variations  in  the  corrected  m^  residuals 
with  event  location  are  indeed  associated  with  variations 
in  the  subsurface  geologic  structure  across  the  site,  then 
it  should  be  possible  to  contour  the  variations  in  station 
corrections  as  a  function  of  source  location.  This  in  fact 
appears  to  be  the  case,  although  the  limited  spatial  distri¬ 
bution  and  precision  of  the  available  explosion  data  makes 
it  difficult  to  do  anything  more  than  outline  the  broad 
trends  at  the  present  time.  For  example,  Figure  25  shows 
the  average  station-corrected  residuals  as  a  function  of 
event  location,  computed  using  the  five  stations  in  the 


51 


'J-  V 


Degrees  North  Latitude 


50.10 


Degrees  East  Longitude 


Figure  25.  Contours  of  mean  residuals  across  the  Shagan 
test  site  derived  from  a  group  of  stations  with 
an  average  recording  azimuth  of  3®.  Contour  is 
in  m^  units.  Solid  circles  are  event  locations 
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northern  azimuth  (i.e.,  340°  <  8  <  20°)  which  recorded  80 
percent  or  more  of  the  selected  explosions.  The  computed 
95  percent  confidence  intervals  about  the  mean  values  shown 
on  this  figure  average  about  t  0.08  magnitude  units,  which 
provides  some  idea  of  the  resolution  of  these  data.  Rough 
zero  level  contours  have  been  drawn  on  this  figure  to  high¬ 
light  the  general  trends.  It  can  be  seen  that,  as  was  noted 
previously  in  conjunction  with  the  discussion  of  Figure  13, 
the  average  station-corrected  residuals  for  events  in  the 
southwest  portion  of  the  test  site  are  small  and  remarkably 
uniform  over  a  broad  area.  In  contrast  to  this,  the  average 
station-corrected  m^  residuals  for  events  in  the  northeast 
and  central  portions  of  the  test  site  are  relatively  large 
and  vary  by  more  than  0.3  units  m^  over  distances  of  less 
than  5  km.  However,  despite  the  fact  that  this  spatial  var¬ 
iation  is  rapid,  it  is  nevertheless  a  fairly  systematic 
function  of  location,  as  can  be  seen  with  reference  to  the 
zero  level  contour  drawn  through  that  area.  For  purposes  of 
completeness,  the  station-corrected  m^  residuals  for  each  of 
the  15  recording  stations  located  in  the  northern  azimuth 
window,  including  those  used  in  deriving  Figure  25,  are 
shown  as  a  function  of  event  location  in  the  central  and 
northeast  portions  of  the  test  site  in  Appendix  A,  where  the 
zero  level  contour  from  Figure  25  has  been  superimposed  on 
each  station  plot  as  a  measure  of  the  internal  consistency 
of  the  complete  data  set. 

The  variation  in  the  average  station-corrected  m^ 
residuals  as  a  function  of  event  location  computed  using 
the  eight  stations  in  the  northwest  azimuth  (i.e.,  290°  <  8  < 
305°)  which  recorded  80  percent  or  more  of  the  selected  ex¬ 
plosions  is  shown  in  Figure  26.  In  this  case,  due  to  the 
larger  sample  size,  the  computed  95  percent  confidence  in¬ 
tervals  associated  with  the  mean  values  shown  on  the  figure 
average  about  ±  0.06  magnitude  units.  It  can  be  seen  from 
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Thus,  for  this  event,  the  apparent  m.  residuals  for  stations 
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It  follows  that,  for  example  if  =  3N2,  then 
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which  illustrates  how  the  actual  bias  in  the  northwest 
azimuth  can  be  partitioned  to  give  a  large  apparent  bias  in 
the  northern  azimuth.  By  the  same  line  of  reasoning,  suppose 
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for  a  second  event  that  m^  *  nfc2  +  Am'  where  again  mj,2  = 
mb2  and  =  3N2.  In  this  case 
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as  it  should.  Thus,  with  reference  to  Figure  13,  it  is  the 
difference  between  the  average  residual  change  in  various 
azimuth  windows  that  is  well  constrained,  and  not  the  abso¬ 
lute  value  of  the  residual  change  in  any  one  azimuth.  That 
is,  the  absolute  Am^  scale  on  that  figure  can  be  shifted  up 
or  down  by  changing  the  azimuthal  distribution  of  the  sta¬ 
tions  in  the  network.  It  follows  that  one  cannot  conclude 
from  a  comparison  of  Figures  25  and  26  that  the  variability 
is  strongest  at  stations  in  the  northern  azimuth.  What 
these  figures  do  demonstrate  is  that  the  variability  with 
source  location  is  quite  regular  and  most  pronounced  between 
the  northeast  and  central  portions  of  the  test  site,  indi¬ 
cating  that  a  change  in  the  subsurface  geology  must  occur  in 
this  vicinity. 

3.2  P  WAVE  TRAVEL-TIME  DATA 

Analyses  of  travel-time  data  recorded  from  explo¬ 
sions  at  the  Shagan  River  test  site  provide  another  potential 
means  for  investigating  test  site  variability .  In  constrast 
to  the  m^  data  discussed  above,  most  of  the  stations  report 
arrival  time  data  for  all  the  selected  explosions  which 
occurred  during  the  period  in  which  they  were  operational. 

In  fact,  on  the  average,  80  of  the  94  stations  discussed  in 
Section  II  report  arrival  times  for  any  given  explosion. 

Our  analysis  of  these  data  was  initiated  by  transforming  the 
reported  arrival  times  into  estimates  of  the  travel  times  by 
subtracting  the  JED  event  origin  times  obtained  by  Marshall 
et  al.  (1984).  Event-station  distances  were  then  computed 
using  the  JED  explosion  epicenters  listed  in  Table  1  and 


travel-time  residuals  were  computed  by  subtracting  the  travel 
times  predicted  by  the  Herrin  1968  P  wave  travel-time  curves 
(Herrin  et  al . ,  1968)  from  the  observed  travel  times.  Then, 
by  analogy  with  the  m^  analysis  described  above,  the  mean 
travel-time  residual  was  determined  for  each  station,  and 
subtracted  from  the  individual  observed  residuals  at  that 
station  to  obtain  the  variation  of  the  residuals  as  a  func¬ 
tion  of  source  location. 

Unlike  the  m^  residual  data,  the  reduced  travel¬ 
time  residuals  were  found  to  be  quite  small  and  to  show  no 
obvious  correlation  with  event  location.  This  is  illustrated 
in  Figure  27  which  shows  the  reduced  travel-time  residuals 
as  a  function  of  event  location  at  ten  different  stations 
representing  a  wide  range  of  azimuths.  It  can  be  seen  that 
these  station-corrected  travel-time  residuals  generally  fall 
in  the  range  of  ±  0.5  seconds  and  appear  to  vary  randomly 
across  the  test  site.  Similarly,  Figure  28  shows  the  azi¬ 
muthal  distribution  of  the  station-corrected  travel-time 
residuals  for  the  same  central  and  northeast  Shagan  explo¬ 
sions  for  which  m^  data  were  shown  in  Figure  14 .  As  in  the 
plot,  station  azimuth  in  Figure  28  is  measured  clockwise 
from  north,  and  the  circle  corresponds  to  a  corrected  travel¬ 
time  residual  of  zero,  with  positive  residuals  plotting  out¬ 
side  the  circle  and  negative  residuals  plotting  inside  the 
circle  according  to  the  time  scale  shown  on  the  figure. 
Comparing  Figures  28  and  14,  it  can  be  seen  that  the  cor¬ 
rected  travel-time  residuals  do  not  show  the  systematic  pat¬ 
tern  of  the  corresponding  corrected  m^  residuals  and  appear 
to  be  randomly  distributed  about  a  mean  of  zero  at  all  azi¬ 
muths,  independent  of  source  location. 

Thus,  the  observed  travel-time  residuals  do  not 
show  any  clear  correlation  with  the  corresponding  m^  resid¬ 
uals.  This  is  not  a  surprising  observation  and  has,  in  fact, 
been  frequently  noted  by  previous  investigators.  Thus,  for 
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Figure  28.  Aiimuthal  distribution  of  corrected  travel  time  residuals  as 
function  of  event  location;  NE  and  central  events.  Circle 
represents  tero  residual;  1 -second  range  is  denoted. 


example,  Minster  et  si .  (1981)  found  that  variations 
across  MTS  of  tha  order  of  0.4  magnitude  units  were  accom¬ 
panied  by  variations  in  the  travel-time  residuals  of  less 
than  0.4  seconds.  Travel-time  variations  of  this  size  are 
difficult  to  resolve  from  the  rather  imprecise  times  re¬ 
ported  to  ISC ,  particularly  in  view  of  the  fact  that  the 
actual  depths  of  the  explosions  at  Shagan  River  are  not 
known.  In  any  case,  it  does  appear  that  these  P  wave  arri¬ 
val  time  data  will  not  be  useful  for  identifying  any  varia¬ 
tions  in  the  subsurface  geology  across  the  Shagan  River 
testing  area. 


IV.  SUMMARY  AMD  PRELIMINARY  CONCLUSIONS 


4.1  SUMMARY 

The  investigation  summarized  in  this  report  has 
centered  on  a  preliminary  analysis  of  seismic  variability  at 
the  Shagan  River  nuclear  test  site.  Specifically,  large 
samples  of  teleseismic  P  wave  amplitude  and  arrival  time 
data  recorded  from  explosions  at  this  test  site  have  been 
collected  and  statistically  analyzed  in  an  attempt  to  define 
any  systematic  trends  which  correlate  with  source  location. 

The  teleseismic  P  wave  data  base  which  has  been 
assembled  for  this  project  was  described  in  Section  II  where 
the  distributions  of  the  data  with  respect  to  source  and 
station  parameters  were  documented  in  detail.  The  data  base 
is  composed  of  individual  readings  and  P  wave  arrival 
times  recorded  by  a  selected  network  of  94  worldwide  receiver 
stations  from  a  sample  of  52  Shagan  River  underground  nuclear 
explosions  which  have  been  assigned  m^  values  of  5.5  or 
greater.  These  seismic  data,  together  with  the  refined  ex¬ 
plosion  epicenters  and  origin  times  determined  by  Marshall 
et  al.  (1984) ,  constitute  the  data  set  which  was  used  to  in¬ 
vestigate  test  site  variability. 

A  systematic  analysis  of  these  data  was  described 
in  Section  III  where  the  LSMF  statistical  analysis  procedure 
was  used  to  derive  network -averaged  event  magnitudes  and 
average  m^  station  correction  factors  for  explosions  at  the 
Shagan  River  test  site.  The  variations  of  these  station 
correction  factors  with  source  location  within  the  test  site 
were  then  carefully  analyzed  and  it  was  shown  that  large 
changes  of  the  order  of  0.5  units  can  occur  between  rela¬ 
tively  closely  spaced  explosions.  Several  alternate  explan¬ 
ations  of  these  m^  variations  were  then  tested  and  critically 
evaluated  through  comparisons  of  data  recorded  at  common 
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stations  from  selectsd  pairs  of  explosions.  This  was  fol¬ 
lowed  by  a  preliminary  analysis  of  the  corresponding  P  wave 
arrival  time  data. 

4.2  PRELIMINARY  CONCLUSIONS 

The  analyses  summarized  above  support  the  following 
preliminary  conclusions  regarding  seismic  variability  at  the 
Shagan  River  nuclear  test  site. 

(1)  Tele seismic  data  provide  strong  evidence 
of  systematic  geophysical  variations  within 
the  Shagan  River  testing  area.  In  particular, 
corrected  single  station  residuals  for 
explosions  in  the  southwest  portion  of  the 
test  site  appear  to  be  random  at  any  given 
azimuth  and  quite  similar  from  event  to  event, 
while  there  are  large  (0.5  units  n^)  varia¬ 
tions  in  the  corrected  nt  residuals  with  azi¬ 
muth  between  explosions  in  close  proximity  in 
the  northeast  and  central  portions  of  the  test 
site . 

(2)  The  observed  azimuthal  variability  in  the 
station  corrections  with  source  location  in¬ 
dicates  that  the  network-averaged  m^  values, 
and  corresponding  yield  estimates,  for  explo¬ 
sions  at  Shagan  River  will  be  dependent  on  the 
specific  azimuthal  distribution  of  the  sta¬ 
tions  used  to  compute  the  averages.  This  is 
an  additional  source  of  uncertainty  to  be 
considered  in  evaluating  yield  estimates 
determined  from  small  network  averages. 

(3)  The  results  of  detailed  comparisons  of  m^ 
residual  data  at  common  stations  from  selected 
pairs  of  explosions  strongly  suggest  that  the 


observed  variations  in  the  azimuthal  patterns 
of  the  m^  residuals  are  related  to  event  loca¬ 
tion  rather  than  tectonic  release  effects. 

(4)  Preliminary  analyses  indicate  that  the  varia¬ 
tions  in  m^  corrections  for  stations  in  a 
given  azimuth  are  systematic  enough  that  they 
can  be  contoured  as  a  function  of  source  loca¬ 
tion.  This  result  supports  the  hypothesis 
that  the  observed  differences  are  associated 
with  changes  in  the  near-source  P  wave  propa¬ 
gation  paths  to  teleseismic  distances  as  a 
function  of  source  location  within  the  test 
site. 

(5)  In  contrast  to  the  m^  residual  data,  the  re¬ 
duced  P  wave  travel-time  residuals  show  no 
obvious  correlation  with  event  location.  This 
is  interpreted  as  an  indication  that  the  arri¬ 
val  time  data  reported  to  ISC  is  not  precise 
enough  to  resolve  the  rather  small  travel-time 
variations  which  might  be  expected  to  accompany 
the  observed  m^  variations.  In  any  case,  it 
does  not  appear  that  the  ISC  travel  time  data 
will  be  useful  for  identifying  variations  in 
the  subsurface  geology  across  the  Shagan  River 
test  site. 
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APPENDIX  A 

Station-corrected  m^  residuals  as  a  function  of 
source  location  for  events  in  central  and  northeastern 
Shagan  recorded  by  all  individual  stations  within  the  azi¬ 
muthal  window  having  an  average  of  3°.  Stations  denoted 
with  (*)  were  used  to  compute  mean  residuals  for  contour 
plot  in  Figure  25.  Shown  here  is  the  zero  m^  residual 
contour  line  from  Figure  25.  Area  displayed  is  49.92  - 
50.10  degrees  North  latitude  by  78.90  -  79.05  degrees  East 
longitude. 
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